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ABSTRACT 


B.  A.  Kimball,  D.  L.  Kittock,  T.  J.  Henneberry,  and  L.  A.  Bariola.     Computer  Pro- 
gram for  Predicting  the  Reduction  of  Pink  Bollworm  Populations  by  Chemically 
Terminating  Late-Season  Bolls  in  Cotton.     U.S.  Department  of  Agriculture, 
Agricultural  Research  Service  ARS  W-48,  48  pp.,  September  1977. 


Diapausing  pink  bollworms,  Pectinophora  gossypiella  (Saunders),  develop  in 
unharvested  immature  cotton  bolls  late  in  the  growing  season.  Consequently, 
eliminating  these  bolls  by  chemical  termination  has  proven  to  be  an  effective 
control  method.     However,  if  the  chemical  treatment  is  applied  too  early  in  the 
season,  some  bolls  that  should  have  contributed  to  the  yield  may  also  be  elimi- 
nated; so  timing  of  the  chemical  application  is  important.     Therefore,  a  computer 
model  was  developed  for  predicting  the  effects  on  pink  bollworm  populations  and 
on  cotton  yields  of  terminating  boll  development  with  chemicals  applied  on  any 
given  date  during  the  growing  season. 

The  program  listing  is  presented  along  with  flow  diagrams  and  a  description 
of  how  it  works.     The  program  requires  field  data  on  the  daily  boll  set  rate  for 
a  particular  cotton  variety,  on  the  rate  at  which  the  bolls  mature  during  the 
season,  on  the  reduction  of  boll  set  by  a  particular  chemical,  and  on  the  sea- 
sonal patterns  of  bollworm  infestation  and  diapause  transition.     An  example  is 
presented  of  the  results  obtained  for  a  specific  case. 
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This  paper  contains  the  results  of  research  only.     Mention  of 
a  pesticide  in  this  paper  does  not  constitute  a  recommendation 
for  use  by  the  U.S.  Department  of  Agriculture  nor  does  it  imply 
registration  under  FIFRA  as  amended. 


Trade  names  and  the  names  of  commercial  companies  are  used 
in  this  publication  solely  to  provide  specific  information. 
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COMPUTER  PROGRAM  FOR  PREDICTING  THE  REDUCTION  OF 
PINK  BOLLWORM  POPULATIONS  BY  CHEMICALLY  TERMINATING 
LATE-SEASON  BOLLS  IN  COTTON 


By  B.  A.  Kimball,  D.  L.  Kittock,  T.  J.  Henneberry, 
and  L.  A.  Bariola-'- 

THEORY  OF  FIELD  METHOD 

The  larvae  of  the  pink  bollworm,  Pectinophora  gossypiella   (Saunders) ,  feed 
in  green  cotton  bolls.     When  the  larvae  mature,   they  "cut  out"  and  drop  to  the 
ground,  where,  in  summer,  they  pupate.     However,  starting  in  September,  progres- 
sively larger  numbers  of  the  larvae  do  not  pupate,  but  instead  enter  the  dia- 
pause or  overwintering  stage.     These  diapausing  larvae  develop  in  late-season 
green  bolls,  which  usually  do  not  mature  by  harvesttime.     These  bolls  can  be 
eliminated  by  early  defoliation,  desiccation,  or  early  termination  of  irrigation. 
However,  these  practices  affect  maturation  of  all  bolls,  and,  consequently,  cot- 
ton yield  and  quality  are  reduced.     A  less  severe  treatment  was  proposed  by 
Kittock  et  al.    (4,  5,  6)^  and  Mauney  et  al.    (7).     They  applied  growth  regulators 
to  eliminate  the  late-season  green  bolls,  but  allowed  bolls  already  set  to  con- 
tinue their  development  and  be  harvested.     Bariola  et  al.    (1)  showed  that  appli- 
cations of  certain  of  these  chemicals  reduced  the  number  of  immature  green  bolls 
and  the  number  of  diapausing  pink  bollworm  larvae  with  little  or  no  reduction  of 
yield.     The  purpose  of  this  computer  program  is  to  predict  the  diapausing  pink 
bollworm  population,  the  reduction  in  cotton  yield  when  such  chemicals  are  ap- 
plied to  cotton,  and  the  optimum  time  to  apply  these  chemicals. 

DESCRIPTION  OF  COMPUTER  MODEL 

Figures  1  and  2   (Kittock,  unpublished)  illustrate  how  various  cotton  vari- 
eties set  bolls  during  a  season.     Figure  3  shows  the  length  of  time  required  for 
the  bolls  to  mature.     The  bolls  mature  faster  at  the  higher  temperatures  present 
early  in  the  season,  and  slower  later  in  the  season.     However,  the  bolls  set  on 
any  particular  day  do  not  all  mature  at  exactly  the  same  rate  but  instead  are 
normally  distributed  about  a  mean  rate  (fig.  4). 

Thus,  the  curve  shows  when  50  percent  of  the  bolls  set  on  any  particular 
day  will  have  matured,  and  a  normal  distribution  is  superimposed  on  this  curve. 
By  summing  the  contributions  from  all  previous  days,  the  number  of  bolls  maturing 
on  a  particular  day  is  computed.     By  accumulating  the  number  of  bolls  maturing 
each  day  through  the  season,  the  final  yield,  in  bolls  per  acre,  is  obtained. 


■•■Kimball  is  a  soil  scientist,  U.S.  Water  Conservation  Laboratory;  Kittock 
is  an  agronomist,  University  of  Arizona  Cotton  Research  Center;  and  Henneberry 
and  Bariola  are  entomologists.  Western  Cotton  Research  Laboratory.     All  authors 
are  with  the  Agricultural  Research  Service,  Phoenix,  Ariz. 

Italic  numbers  in  parentheses  refer  to  Literature  Cited,  p.  17. 
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Figure  1. — Boll  set  versus  time  during  season  for  DPL  16  cotton.     The  curve  is 
based  on  1973  data  from  Parker,  Ariz.     The  arrows  indicate  the  knot  posi- 
tions, which  delimit  the  subranges  used  with  the  cubic  spline  smoothing. 
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Figure  2. — Boll  set  versus  time  during  season  for  AZ  64  cotton.     The  curve  is 

based  on  1975  data  from  Phoenix,  Ariz. 
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Figure  4. — Distribution  curves  used  to  model  the  differing 
physiological  maturation  rates  of  cotton  bolls. 
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Figure  5. — Reduction  in  boll  set  of 
cotton  caused  by  a  treatment  with 
2,4-D  plus  chlorf lurenol  versus 
time  from  day  of  treatment. 
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Figure  6. — Reduction  in  boll  set  of 
cotton  caused  by  a  treatment  with 
TD-1123  plus  chlorf lurenol  versus 
time  from  day  of  treatment. 


Figures  5  and  6   (Kittock,  unpublished)  show  how  boll  set  is  affected  as  a 
function  of  days  from  time  of  chemical  treatment. 

A  new  boll  set  pattern  is  obtained,  but  the  yield  is  obtained  by  the  same 
procedure  as  explained  above  for  no  chemical  treatment.     The  chemical  names  of 
the  compounds  used  are:     2,4-D,   (2 , 4-dichlorophenoxy)  acetic  acid;  chlorf lurenol , 
methyl  2-chloro-9-hydroxyf luorene-9-carboxylate ;  and  Pennwalt  TD-1123,  potassium 
3, 4-dichloroisothiazole  5-carboxylate. 
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The  number  of  pink  bollworm  larvae  is  more  difficult  to  compute  but  can  be 
obtained  using  infestation  data  and  some  assumptions.     Figures  7  and  8  (Vail  et 
al. ,  8)  show  the  typical  number  of  larvae  per  green  boll  as  a  function  of  time 
during  the  season. 

Thus,  if  we  can  compute  the  number  of  green  bolls  in  a  field,  for  similar 
infestation  conditions  we  can  compute  the  number  of  larvae  in  the  field  as  a 
function  of  time  during  the  season.     Figure  9  shows  the  fraction  of  the  larvae 
going  into  diapause  as  a  function  of  the  time  of  year  (Bariola,  unpublished). 
Thus,  once  the  number  of  larvae  in  a  field  is  obtained,  the  number  of  diapaused 
larvae  can  also  be  obtained. 

First,  the  number  of  green  bolls  is  computed.     Only  green  bolls  are  suscep- 
tible to  the  larvae.     The  age  span  of  the  susceptible  bolls  was  found  to  be  10 
to  28  days  early  in  the  season.     Since  bolls  mature  more  slowly  later  in  the 
season,  we  assumed  that  late-season  bolls  will  require  more  than  10  days  to  be- 
come susceptible  and  that  they  will  remain  susceptible  beyond  the  age  of  28  days. 
Therefore,  the  first  age  of  susceptibility  is  computed  as  lYNG^  =  10.     *  (AVDA/ 
45.)  and  the  final  age  as  lOLD  =  28.  *  (AVDA/A5.)  where  AVDA  is  the  average  num- 
ber of  days  required  for  the  particular  day's  bolls  to  mature.     The  normalizing 
value  of  45.  comes  from  the  fact  that  about  45  days  are  required  for  early  sea- 
son bolls  to  mature  (fig.  3). 

Just  as  there  is  a  normal  distribution  of  maturation  rates  that  must  be 
used  to  compute  the  lengths  of  time  required  to  reach  full  maturation,  it  is 
plausible  that  there  is  also  an  age  distribution  for  becoming  susceptible  to 
pink  bollworm  and  another  for  ceasing  to  be  susceptible  to  pink  bollworm.  As 
shown  in  figure  4,  for  maturation  at  a  mean  of  50  days,  90  percent  of  the  bolls 
mature  between  41  and  59  days,  which  corresponds  to  a  standard  deviation  of  5.47 
days.     This  curve  is  based  on    field  data  (Kittock,  unpublished).     The  bolls 
with  this  maturation  distribution  were  all  set  on  the  same  day,  so  an  initial 
rectangular  "pulse"  has  become  smeared  into  a  flattened  bell  after  a  mean  of  50 
days.     It  is  plausible  that  the  standard  deviation  of  this  boll  distribution 
changes  in  proportion  to  the  mean.     Thus,  we  compute  the  standard  deviation, 
SIGMA,  for  maturation  distribution  as  SIGMA  =  (AVDA/50.)  *5.47.     Similarly,  we 
compute  its  standard  deviation  for  the  first  susceptibility  to  pink  bollworm 
attack  as  SIGMA  =  (IYNG/50.)  *5.47,  and  for  cessation  of  susceptibility  as 
SIGM  =  (IOLD/50.)  *5.47.     The  resultant  distributions  for  means  of  10  and  28 
days  are  plotted  in  figure  4.     The  10-day  curve  shows  98  percent  become  suscep- 
tible between  8  and  12  days  of  age,  and  the  28-day  curve  shows  93  percent  cease 
to  be  susceptible  between  23  and  33  days  of  age. 

From  the  rate  of  maturation  and  its  distribution,  as  described  in  the  pre- 
vious two  paragraphs,  the  number  of  susceptible  green  bolls  can  be  computed. 
This  is  done  by  establishing  three  computational  loops.     The  outermost  loop  is 
for  days  through  the  season  indexed  by  the  symbol  IDATE.     The  next  loop  sums  the 
contribution  to  the  green  boll  population  present  on  IDATE  from  all  days  previ- 
ous to  IDATE.     The  second  loop  is  indexed  by  JDATE.     The  innermost  loop  indexed 
by  IGRN  computes  the  number  of  IGRN-days-old  bolls  that  were  set  on  JDATE  and 
are  in  the  susceptible  boll  population  present  on  IDATE. 


^The  symbolic  names  of  variables  used  in  the  FORTRAN  computer  program  are 
listed  and  defined  in  the  Nomenclature  section  of  the  Appendix,  p.  18. 
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Figure  7. — Curve  showing  the  increasing  infestation  of  green  cotton  bolls  with 
pink  bollworm  larvae  during  the  growing  season  for  a  field  that  received 
no  application  of  insecticides. 
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Figure  8. — Curve  showing  the  increasing  infestation  of  green  cotton  bolls  with 
pink  bollworm  larvae  during  the  growing  season  for  a  field  that  received 
applications  of  insecticides. 
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Figure  9. — Curve  shoving  the  percentage  of  pink  bollworm  larvae  going  into 
diapause  versus  time  during  the  growing  season  when  the  larvae  infested 
a  green  cotton  boll. 

As  shown  in  flow  diagrams  in  the  Appendix,  the  number  of  lYNG  bolls  is  com- 
puted first.     The  standard  deviation  of  this  group  of  bolls  is  obtained  from 
SIGM  =  (IGRN/50.)  *5.45  as  explained  above,  where  IGRN  =  lYKG  the  first  time 
through  the  loop.     Then  utilizing  an  error  function  table,  the  fraction  (F) , 
which  these  bolls  represent  of  those  set  on  JDATE,  is  computed  as  the  area  under 
the  normal  distribution  curve  for  this  particular  mean  boll  age.     The  integra- 
tion limits  go  from  JDIF  -  0.5  to  JDIF  +0.5  days,  where  JDIF  =  JDATE  +  IGRN  - 
IDATE,  the  difference  between  IDATE  and  the  mean  age  of  these  bolls.     As  the  in- 
nermost loop  is  repeated  for  each  green  boll  age  (IGRN) ,  F  is  obtained  for  each 
age  and  is  miultipled  by  the  boll  set  num.ber  for  JDATE  to  get  the  number  of  each 
age  of  bolls  contributed  by  JDATE.     The  same  F's  are  used  to  obtain  the  green 
boll  numbers  for  when  chemical  treatments  were  applied  by  using  the  corresponding 
reduced  boll  set  number  for  JDATE  obtained  earlier  to  account  for  chemical 
treatments.     The  green  boll  numbers  without  and  with  chemical  treatments  are 
stored  in  arrays,  which  are  indexed  by  IGEi^.     Each  element  stores  the  accumula- 
tion of  IGRN  day-old  bolls  from  previous  days   (JDATE)  to  the  population  present 
on  IDATE. 

After  the  population  of  green  bolls  has  been  obtained  for  IDATE,  the  amount 
of  new  infestation  of  pink  bollworms  is  computed  for  the  case  of  no  chemical 
treatment.     Using  the  infestation  curve  (figs.   7  and  8),  the  number  of  new  pink 
bollworm  larvae  is  computed  as  the  sum  of  the  average  fraction  for  IDATE  (di- 
rectly from  the  infestation  curve) ,  times  the  number  of  green  bolls  that  have 
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become  newly  susceptible,  plus  the  change  in  infestation  since  the  previous  day, 
times  the  number  of  all  other  ages  of  susceptible  green  bolls.     This  procedure 
is  valid  only  for  the  check  field  with  no  chemical  termination  treatment  because 
the  infestation  curves  came  from  such  untreated  fields,  which  have  a  different 
green  boll  age  distribution  than  the  treated  fields. 

Once  the  total  number  of  new  pink  bollworm  larvae  (PNEW)  has  been  computed 
for  the  case  of  no  chemical  termination  treatment,  an  adjustment  is  made  to  ob- 
tain the  number  of  new  pink  bollworm  larvae  (CPNEW)  for  the  case  where  there  is 
a  treatment.     First,  we  assume  there  is  an  equal  probability  for  all  susceptible 
green  bolls  to  be  infested  with  new  larvae.     The  amount  of  new  infestation  per 
boll  is  PNEW/TGBLS,  where  TGBLS  is  the  total  number  of  green  bolls.     The  number 
of  new  larvae  in  a  particular  age  of  boll  (PVBVD)  is  then  simply  the  number  of 
bolls  of  that  age  times  PNEW/TGBLS.    We  next  assume  that  when  the  chemical 
treatment  removes  green  bolls,  the  pink  bollworm  population  is  reduced  propor- 
tionately and  infestation  of  remaining  bolls  does  not  increase.     Kittock  et  al. 
(4)  have  shown  this  is  a  reasonable  assumption.     Thus,  assuming  that  the  pink 
bollworm  population  is  reduced  proportional  to  the  green  boll  reduction,  CPNEW 
is  computed  as  the  sum  of  over  boll  age  of  PVBVD  times  the  ratio  of  number  of 
green  bolls  with  to  without  chemical  treatment  of  each  boll  age. 

Finally,  after  the  amounts  of  new  infestation  without  and  with  a  chemical 
termination  treatment  have  been  computed  for  IDATE,  the  fraction  going  into  a 
diapause  is  computed  from  the  diapause  curve  (fig.  9).     Then,  the  total  number 
going  into  diapause  is  accumulated  through  the  season  as  IDATE  steps  along  day 
by  day. 


INPUT  CURVES  AND  FIGURES  FROM  FIELD  DATA 

Since  various  types  of  cotton  have  various  patterns  of  setting  bolls  through 
a  season,  and  since  different  chemical  treatments  affect  boll  set  differently, 
we  decided  at  the  outset  that  the  various  input  curves  should  be  represented  by 
individual  subroutines.     This  method  allows  the  same  program  to  be  used  while 
only  substituting  one  subroutine  for  another  as  one  changes  from  one  case  to 
another. 

The  subroutines  are  all  similar  in  that  all  have  two  arguments.     The  first 
is  a  number  or  a  date  given  to  the  subroutine  by  the  main  calling  program,  and 
the  second  is  a  value  computed  by  the  subroutine  and  returned  to  the  main  program. 

The  data  in  all  figures  (except  fig.  4)  have  been  approximated  using  smooth 
curves.     The  method  of  smoothing  was  that  of  cubic  splines   (Kimball,  2,   5),  which 
fits  cubic  polynomials  to  subranges  of  the  data*.     The  end  points  of  the  subranges 
are  called  knots  and  are  indicated  by  arrows  in  the  figures.     Each  subrange  of 
data  has  a  cubic  equation  that  was  fitted  subject  to  the  constraints  that  it  must 
have  the  same  function  value  and  the  same  first  and  second  derivative  at  the  knot 
where  it  meets  the  cubic  equation  for  the  adjacent  subrange.     The  values  of  the 
coefficients  for  a  subrange  of  a  particular  curve  can  be  obtained  by  looking  at 
the  subroutine    for  the  particular  curve  in  the  program  section  of  the  Appendix, 
page  20.     First,  sort  through  the  IF  statements  to  find  the  particular  subrange. 
Then,  the  numbers  for  the  CUBIC  function  are  the  constant,  linear,  quadratic,  and 
cubic  coefficients,  respectively.     Often,  the  first  and  last  subrange  were  delib- 
erately set  to  an  appropriate  constant  such  as  zero.     Other  smoothing  methods  than 


8 


cubic  splines  can  also  be  used,  or  even  tabular  data,  so  long  as  the  subroutines 
have  the  correct  name  and  argument  fomat. 

The  curves  presented  here  were  used  to  represent  types  of  data  rather  than 
to  represent  particular  sets  of  data  as  closely  as  possible.     Therefore,  only 
the  curves  are  presented,  and  not  the  original  data  points  upon  which  they  are 
based.     Thus,  the  computed  results  can  be  regarded  as  typical  of  what  to  expect 
when  the  input  curves  are  of  the  type  shown. 

Boll  Set 

Figure  1  shows  the  boll  set  pattern  for  DPL  16  (Deltapine  16)  cotton  ob- 
tained in  1973  at  Parker,  Ariz.     The  August  slump  in  boll  set  followed  by  the 
setting  of  late-season  bolls,  many  of  which  do  not  mature,  makes  this  variety 
of  cotton  a  particularly  good  candidate  for  reducing  pink  bollworm  populations 
by  chemical  termination  of  the  late-season  bolls. 

Figure  2  shows  the  boll  set  pattern  for  AZ  64  cotton  based  on  1975  data 
from  Phoenix.     This  pattern  represents  an  opposite  type  of  cotton  and  growing 
conditions  from  the  DPL  16  of  figure  1,  because  a  large  portion  of  the  bolls  are 
set  in  late  August. 

SUBROUTINE  BOLSET  (IDATE,  B)  is  used  to  compute  the  bollset  where  IDATE  is 
the  Julian  date  of  flowering  given  by  the  main  calling  program,  and  B  is  the 
boll  set  per  acre  per  day  computed  by  BOLSET  and  returned  to  the  main  program. 

Doll  Maturation  Time 

Figure  3  shows  the  length  of  time  for  maturation  required  for  50  percent  of 
the  bolls  set  at  the  time  during  the  season  shown  on  the  axis.     This  curve  is 
based  on  1967  data  from  Phoenix. 

SUBROUTINE  BOLMAT  (IDATE,  AVDA)  is  used  to  compute  the  length  of  time  for 
boll  maturation,  where  IDATE  is  the  Julian  date  of  flowering  given  by  the  calling 
program,  and  AVDA  is  the  number  of  days  required  for  50  percent  of  the  bolls  set 
on  IDATE  to  mature.     AVDA  is  computed  by  BOLMAT  and  returned  to  the  calling 
program. 

Distribution  of  Boll  Maturation  Rates 

Figure  4  shows  three  maturation  rate  distribution  curves  used  by  the  pro- 
gram.    These  curves  are  not  input  to  the  program,  but  instead  are  computed  by 
the  program  from  values  of  the  mean  days  and    standard  deviations.     The  50-day 
curve  is  based  on  actual  data  from  Phoenix  in  1967.     The  data  shows  that  5  per- 
cent of  the  bolls  matured  in  41  days;  50  percent  in  50  days;  and  95  percent,  in 
59  days.     The  distribution  appeared  normal  with  a  mean  of  50  days  and  a  standard 
deviation  of  5.47  days.     The  other  distributions  are  computed  by  making  the 
standard  deviation  proportional  to  the  mean  age  of  the  population,  that  is, 
a  (y/50. )*5.47.     If  a  different  maturation  rate  distribution  is  to  be  used,  the 
main  program  will  have  to  be  changed,  but  this  can  be  done. 
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Chemical  Termination  Treatments 


Figures  5  and  6  show  the  effects  of  2,4-D  plus  chlorf lurenol  and  of  TD-1123 
plus  chlorf lurenol ,  respectively,  on  the  boll  set  of  cotton.     Both  treatments 
affect  bolls  set  up  to  9  days  prior  to  the  day  of  treatment.     TD-1123  stops  boll 
set  almost  immediately,  whereas  2,4-D  allows  some  bolls  to  set  even  20  days  after 
treatment . 

SUBROUTINE  CHEM  (JDIF,  FRACT)  is  used  to  compute  the  effect  of  the  chemical 
treatment  on  boll  set,  where  JDIF  is  the  number  of  days  from  the  date  of  chemical 
treatment  given  by  the  calling  program,  and  FRACT  is  the  fraction  of  bolls  re- 
maining.    FRACT  is  computed  by  CHEM  and  returned  to  the  calling  program. 

Infestation 

Figures  7  and  8  show  the  buildup  of  pink  bollworm  larvae  through  the  grow- 
ing season  for  fields  that  had  not  received  and  had  received,  respectively,  ap- 
plications of  insecticides  for  pink  bollworm  control.     The  numbers  of  larvae 
are  much  higher  during  August  and  September  in  the  no  insecticide  fields.  The 
last  insecticide  application  was  made  on  17  September  on  the  treated  field. 

SUBROUTINE  INFEST  (IDATE,  FNST)  is  used  to  compute  the  infestation,  where 
IDATE  is  the  Julian  date  during  the  growing  season  given  by  the  calling  program, 
and  FNFST  is  the  number  of  pink  bollworm  larvae  per  green  boll.     FNFST  is  com- 
puted by  INFEST  and  returned  to  the  calling  program. 

Diapause 

Figure  9  shows  the  percentage  of  pink  bollworm  larvae  going  into  diapause 
versus  the  Julian  date  of  infestation.     The  big  change  occurs  during  September, 
going  from  about  0  to  90  percent. 

SUBROUTINE  DIPAS  (IDATE,  FDIPAS)  is  used  to  compute  the  fraction  going  into 
diapause,  where  IDATE  is  the  Julian  date  of  infestation  given  by  calling  program, 
and  FDIPAS  is  the  fraction  of  pink  bollworm  larvae  going  into  diapause  that  in- 
fested on  IDATE.     FDIPAS  is  computed  by  DIPAS  and  returned  to  the  calling  program. 

INSTRUCTIONS  FOR  RUNNING  PROGRAM  FROM  COMPUTER  CONSOLE 

The  FORTRAN  program  was  written  for  use  on  a  Data  General  Corporation  Nova 
800  Jumbo  Minicomputer  with  32K  of  memory,  movable  head  disk  storage,  and  a  line 
printer.     First,  the  main  program  and  all  functions  and  subroutines  are  trans- 
ferred to  disk  files  and  compiled.     Then  the  main  program  and  functions,  and 
those  subroutines  needed  for  a  particular  run  are  loaded  and  executed.     The  pro- 
gram will  respond  with  a  question  about  output  files.     If  the  user  responds  with 
a  12,  the  line  printer  will  be  used  for  output.     If  the  user  responds  with  a  2, 
the  output  will  be  stored  in  the  COTRSLT  file  by  the  program.     Next,     the  com- 
puter will  ask  the  number  of  runs.     The  answer  should  be  the  number  of  chemical 
treatment  dates  with  which  this  particular  subroutine  package  is  to  be  run. 
The  computer  will  then  request  chemical  treatment  dates,  one  for  each  run,  in  a 
day,  month  format.     Finally,  it  will  compute  the  numbers  of  cotton  bolls  and 
pink  bollworms  for  each  day  of  the  season  and  output  the  results. 
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EXAMPLE  OF  PRINTED  RESULTS 


Figure  10  (p.   12)  shows  an  example  of  the  printed  results  obtained  from  the 
program  for  a  run  with  DPL  16  cotton  (fig.  1)  that  received  no  insecticide  (fig. 
7)  but  was  treated  with  TD-1123  plus  chlorf lurenol  (fig.   6)  on  Julian  date  244, 
as  indicated  by  the  top  two  heading  lines.     The  next  line  gives  the  date  and  time 
of  the  computer  run.     DA,  MO,  and  DATE  are  the  day,  month,  and  Julian  date 
through  the  growing  season.     MATR  DATE  is  the  day  from  fig.   3  when  50  percent 
of  the  bolls  set  on  the  Julian  date  will  mature.     IFS  is  the  number  of  pink  boll- 
worm  larvae  infesting  100  green  bolls  (fig.  7),  and  DIPS  is  the  percentage  of 
newly  infesting  larvae  going  into  diapause   (fig.   9).     BOLL  SET  is  the  number  of 
bolls  set  per  acre  per  day  (fig.   1)  for  an  untreated  field  (CK)  and  for  a  field 
treated  with  TD-1123  plus  chlorf lurenol  (TMT) .     GREEN  BOLLS,  HARVEST,  AND  DIA- 
PAUSE PBW  are  the  number  of  susceptible  green  bolls,  the  number  of  mature  har- 
vestable  bolls,  and  the  number  of  pink  bollworm  larvae  in  diapause,  respectively, 
per  0.01  acre  for  untreated  (CK)  and  treated  (TMT)  fields. 

Figure  11  (p.   16)     is  a  graph  of  the  numerical  data  shown  in  figure  10.  The 
infestation  (IFS)  and  diapause  (DIPS)   curves  came  from  figures  7  and  9.     The  un- 
treated (CK)  boll  set  curve  is  the  same  as  figure  1,  whereas  the  treated  curve 
shows  the  elimination  of  bolls  set  after  about  23  August  for  a  1  September  appli- 
cation of  TD-1123  plus  chlorf lurenol.     The  number  of  susceptible  green  bolls  in 
the  treated  field  declined  relative  to  the  untreated  field  about  10  days  later, 
and  went  to  zero  by  about  the  end  of  September. 

The  results  of  applying  TD-1123  plus  chlorf lurenol  can  be  seen  by  comparing 
the  number  of  diapausing  pink  bollworms  for  the  untreated  (CK)  and  treated  (TMT) 
fields  (figs.  10  and  11).     For  a  harvest  date  of  30  November,  the  untreated  field 
had  1,623  per  0.01  acre,  vzhereas  the  treated  field  had  only  112.     The  yield  pen- 
alty caused  by  not  allowing  late-season  bolls  to  develop  can  be  seen  by  compar- 
ing the  harvest  bolls  for  the  untreated  (CK)  and  treated  (TMT)  fields.     For  the 
30  November  harvest  date,  the  untreated  field  had  4,236  mature  bolls  per  0.01 
acre,  whereas  the  treated  field  had  3,631.     This  amounts  to  a  14.3-percent  loss 
in  yield  but  a  gain  of  93.1  percent  in  the  reduction  in  the  overwintering  pink 
bollworm  population,  as  stated  at  the  end  of  figure  10. 
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Figure  10. — Computer  output  data. 
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Figure  11. — Plot  of  the  computer  output  data  from  figure  10. 
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APPENDIX 


Nomenclature 


Symbol 


Description 


A  -  Constant  coefficient  of  cubic  polynomial. 

AVDA  -  Days  from  flowering  to  when  50  percent  of  bolls  mature. 

B  -  Bolls  set  per  acre  on  IDATE. 

B  -  Linear  coefficient  of  cubic  polynomial. 

BLSET(     )      -  Array  of  bolls  set  per  day  per  acre.     Each  element  is  the  bolls 

set  on  the  day  given  by  the  subscript. 
BOLMAT  -  Subroutine  to  compute  the  number  of  days  from  flowering  to  when 

50  percent  of  the  bolls  mature. 
BOLSET  -  Subroutine  to  compute  the  bolls  set  per  acre  per  day. 

C  -  Quadratic  coefficient  of  cubic  polynomial. 

CBLSET(     )     -  Array  of  bolls  set  per  day  per  acre  when  there  was  a  chemical 
treatment . 

CDPBW  -  Number  of  pink  bollworm  larvae  (per  acre)  in  diapause  from  field 

that  received  chemical  treatment. 

CGBLS(     )      -  Same  as  GBLS  but  for  a  field  with  a  chemical  treatment. 

CHEM  -  Subroutine  to  compute  the  fractional  reduction  in  boll  set  due  to 

a  chemical  treatment  as  a  function  of  days  from  day  of  application. 

CHRVST  -  Harvested  yield  of  mature  cotton  from  field  that  received  a  chem- 

ical treatment. 

COTRSLT  -  Disk  file  for  program  output. 

CPNEW  -  Same  as  PNEW  but  for  a  field  that  received  a  chemical  treatment. 

CUBIC  -  Function  to  evaluate  a  cubic  polynomial. 

CYIELD  -  Number  of  bolls  per  acre  that  mature  on  IDATE  for  a  field  that 

received  a  chemical  treatment. 
D  -  Cubic  coefficient  of  cubic  polynomial. 

DATE  -  Subroutine  to  read  real  time  clock  and  get  date  program  was  run. 

DIES  -  Subroutine  to  compute  fraction  of  pink  bollworms  from  the  infesta- 

tion newly  occurring  on  IDATE  that  will  go  into  diapause. 
DPBW  -  Number  of  pink  bollworm  larvae  (per  acre)  in  diapause. 

E(     )  -  Table  of  error  function  values. 

ERFUN  (Z)       -  Function  to  calculate  error  function. 

EZM  -  Two  times  the  area  under  the  normal  distribution  curve  integrated 

from  0  to  ZM. 

EZP  -  Two  times  the  area  under  the  normal  distribution  curve  integrated 

from  0  to  ZP. 

F  -  Fraction  of  bolls  started  on  JDATE  maturing  on  IDATE.     In  a  later 

portion  of  the  program,  the  fraction  of  bolls  started  on  JDATE  are 
IGRN  days  old  on  IDATE. 

FDIF  -  Change  of  infestation  from  FYEST  on  the  previous  day  to  FNFST  on 

IDATE. 

FDIPAS  -  Fraction  of  IDATE' s  infestation  going  into  diapause. 

FGTIM  -  Subroutine  to  read  real  time  clock  and  get  time  program  was  run. 

FNFST  -  Ratio  of  larvae  to  bolls  on  IDATE. 

FRACT  -  Fraction  of  bolls  surviving  chemical  treatment,  which  were  set  JDIF 

days  from  a  chemical  treatment. 
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FYEST  -  Ratio  of  larvae  to  green  susceptible  bolls  for  the  day  prior  to 

the  day  currently  being  considered. 
GBLS (     )        -  Array  of  green  bolls  present  on  IDATE.     Each  element  contains  the 

number  of  bolls  per  acre  of  the  age  in  days  given  by  the  subscript. 
KRDUCT  -  Percent  reduction  in  harvest  due  to  chemical  treatment. 

HRVST  -  Harvested  yield  of  mature  cotton  in  bolls  per  acre. 

ID  -  Day  of  month.     After  a  call  to  SR  DIES  with  MO  equal  to  month,  ID 

is  the  Julian  date.     If  ID  is  set  equal  to  the  Julian  date,  and 

MO  to  0,  and  a  call  is  made  to  SR  DIES,   ID  becomes  day  of  month 

again. 

IDATE  -  Julian  date.     Outer  DO  loop  index. 

IGRN  -  Index  for  green  boll  age  loops.     IGRN  is  the  green  boll  age  in 

days . 

IHR  -  Hour  program  was  run. 

IM  -  Integerized  Xl-I. 

IMIN  -  Minute  program  was  run. 

INFEST  -  Subroutine  to  compute  the  ratio  of  larvae  to  green  bolls  for  a 

field  that  received  no  chemical  treatment. 
lOLD  -  Oldest  age  at  which  green  bolls  cease  to  be  susceptible  to  pink 

bollworm. 

TOUT  -  Code  for  output;  12  for  line  printer,  2  for  disk  file  called 

COTRSLT. 
IP  -  Integerized  XP. 

IPl  -  lYNG  +  1. 

IRUN  -  Run  number  index  for  run  loop. 

ISEC  -  Second  program  was  run. 

ITITLE(     )     -  Array  of  text  information  to  identify  particular  program  run. 
lYKG  -  Youngest  age  at  which  green  bolls  become  susceptible  to  pink 

bollworm. 

lYYEST  -  Youngest  susceptible  age,   in  days,  of  green  bolls  for  the  day 

prior  to  the  day  currently  being  considered. 
10-110  -  Scaled  and  integerized  output  variables. 

JDATE  -  Inner  loop  Julian  date.     Inner  loop  is  used  to  compute  the  contri- 

bution of  all  the  days  prior  to  IDATE  to  the  number  of  mature  and 
green  bolls  present  on  IDATE. 

JDIF  -  Difference  between  IDATE  and  day  of  chemical  treatment  application. 

Later  in  program,  it  is  the  difference  between  IDATE  and  maturity 
date,  and,  later  yet,  the  difference  between  IDATE  and  suscepti- 
bility dates. 

KCHEM  -  Date  of  chemical  treatment. 

KCSTOR(  )  -  Array  of  Julian  dates  that  the  chemical  treatment  is  to  be  applied. 
KFLWR  -  Earliest  day  in  year  to  expect  flowering  and,  therefore,  the  date 

to  start  IDATE  loop. 

KHRVST  -  Date  of  final  harvest  or  termination  date  for  running  the  program. 

LPT  -  Line  printer. 

MO  -  Month  of  year  or  when  set  equal  to  0  is  code  for  SR  DIES  to  com- 

pute day  and  month  from  Julian  date. 
NAVDA  -  Integerized  AVDA. 

NDATE(3)        -  Array  having  the  day,  month,  year  program  was  run. 

NRUNS  -  Number  of  runs  with  the  sam.e  set  of  subroutines  and  various  chem- 

ical treatment  dates. 

PORTN  -  Ratio  of  new  infestation  with  no  chemical  treatment  to  the  total 

number  of  green  bolls  for  no  chemical  treatment. 
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PNEW  -  Amount  of  new  infestation,  in  larvae  per  acre,  occurring  on  IDATE 

for  a  field  with  no  chemical  treatment. 
PRDUCT  -  Percent  reduction  in  diapaused  pink  bollworrns  due  to  chemical 

treatment . 

PVBVD  -  Number  of  pink  bollworrns  per  boll  per  IGRN  day's  age. 

SIGMA  -  Standard  deviation  for  normal  distribution  of  maturation  rate. 

TCGBLS  -  Same  as  TGBLS  but  for  a  field  that  received  a  chemical  treatment. 

TGBLS  -  Total  number  of  green  bolls  per  acre  of  all  ages  present  in  field 

on  IDATE. 

XM  -  Lower  integration  limit  equivalent  of  XP. 

XP  -  Scaled  ZP  so  that  the  integer  portion  of  XP  will  be  the  subscript 

for  finding  the  correct  error  function  value  in  E(  )  that  corre- 
sponds to  ZP. 

YIELD  -  Number  of  bolls  per  acre  that  mature  on  the  Julian  date  given  by 

IDATE . 

Z  -  Scaled  variable  for  error  function  argument. 

ZM  -  Similar  to  ZP,  but  represents  the  lower  integration  limits  or  the 

beginning  of  the  JDATE  day. 
ZP  -  Normalized  Z  variable  for  axis  of  normal  distribution  curve.  It 

is  the  upper  integration  limit  representing  the  end  of  the  JDATE 

day. 


Flow  Diagrams 


Figure  12  (pp.  21-32)  consists  of  several  pages  of  flow  diagrams  for  the 
main  program,  each  page  being  a  small  portion  of  a  larger  picture.     To  aid  in 
seeing  the  larger  picture  of  computational  sequence  for  the  whole  program,  pho- 
tocopy the  pages  and  lay  them  end  to  end  so  that  the  circles  with  1,  2,  3,  4,  5, 
7,  8,  12,  13,  14,  and  15  coincide.     Then  draw  lines  to  join  the  other  circles 
with  the  same  number.     The  flow  diagram  follows  the  program  almost  exactly  line 
by  line. 


FORTRAN  Program 


The  computer  listing  of  the  cotton  yield-pink  bollworm  model  is  given  in 
figure  13  (pp.  33-42).     The  program  is  written  in  Data  General  Corporation 
FORTRAN,  but  the  changes  needed  to  run  on  other  computers  are  small.     One  such 
change  that  probably  will  be  necessary  is  to  eliminate  the  five  lines  that  get 
and  write  the  date  and  time  of  the  computer  run.     Another  possible  change  is 
needed  for  computer  systems  that  do  not  have  the  ACCEPT  and  TYPE  statements  for 
inputting  and  outputting  data  from  a  keyboard  terminal.     For  these  systems,  the 
ACCEPT  statement  will  need  to  be  changed  to  the  appropriate  READ  statements,  and 
the  TYPE  statements  can  be  eliminated.     The  first  ACCEPT  statement  requests  a 
value  for  TOUT,  the  output  device  code.     Other  users  should  set  the  value  of 
lOUT  to  the  code  for  their  own  line  printer,  disk  file,  or  other  output  device. 
Other  deviations  from  standard  FORTRAN,  which  require  changing,  are  unlikely. 
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Figure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. 


START 


INPUT  CODE  FOR 
OUTPUT 
DEVICE  lOUT 


N 


Yf:s 


OPEN   DISK    FILE  FOR 
OUTPUT-  COTRSLT 


INPUT-RUN  INFO-ITITLE 
—  NUMBER  RUNS-NRUNS 


IRUN  =  1 


INPUT  DAY  ,  MONTH  OF 
CHEMICAL  TREATMENT 


CONVERT  TO  JULIAN  DATE 
SR  DIES 

N 

STORE  IN    KCSTOR (  ) 

\ 

I  RUN  =  IRUN  +1 
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Figure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


0 


I  RUN   >  N  RUN  P 


NO 


INITIALIZE 

FIRST 

FLOWERING  DATE, 

KFLWR 

INITIALIZE 

FINAL 

HARVEST  DATE.KHRVST 

I  RUN  =  1 


GET  CHEMICAL  TREATMENT   DATE , 
KCHEM    FROM  STORAGE 


WRITE    RUN  INFO 


GET  TIME    FROM  COMPUTER 
CLOCK    AND  OUTPUT 


WRITE  HEADINGS 


COMPUTE     ERROR    FUNCTION  TABLE 


 ^  

INITIALIZE  HARVEST 
W/0     CHEM  TMT  HRVST    =  0. 

W/       CHEM  TMT       CHRVST   =  0. 




INITIALIZE  NO.  PINK    BW   IN  DIAPAUSE 

W/0    CHEM  TMT  DPBW  =  0. 

W/      CHEM  TMT  CDPBW=0. 
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Figure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


INITIALIZE    yesterday's  YOUNGEST 
SUSCEPTIBLE  AGE 
I Y  YEST  =9 


INITIALIZE    YESTERDAYS  AVERAGE 
FRACTIONAL  INFESTATION 
F  YEST  =  0. 


INITIALIZE  BOLLSET  ARRAYS 
W/O  CHEM  TMT  BLSET{l)  =  0. 
W/     CHEM   TMT     CBLSETII)  =  0. 


IDATE  =  KFLWR 


GET    NUMBER  OF   BOLLS, B, SET  ON 
IDATE  IF  NO  CHEM    TMT.  CALL 
BOLSET  tl DATE  ,  B  ) 


BL  SET  (IDATE 

)  =  B 

JDIF  =  I  DATE  - 

KCHEM 

 '±  

GET    FRACTIONAL   REDUCTION   IN  BOLL 
SET   DUE  TO  A  CHEMICAL  TREATMENT 
CALL    CHEM  {  JDIF,  FR ACT  ) 


CB  L  SET  (I  DATE  ] 

=    B  *  FRACT 
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Figure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


® 


\ 

I  DATE  =  I  DATE  +  1 

\ 

INITIALIZE    YIELDS   OBTAINED  ON 
I  DATE  YIELD  =  0. 

C  Y  I  E  LD  =  0. 


INITIALIZE    ARRAYS    OF  GREEN 
BOLLS    PRESENT  ON  IDATE 
GBLS  (I)  =  0. 
CGBLS  (I  )  =  0. 


J  DATE 


KFLWR 


GET  NO.  OF  DAYS   50%  OF  THE  BOLLS 
SET  ON  J  DATE  WILL  TAKE  TO  MATURE 
CALL    BOLMAT  (JDATE  ,  AVDA) 


CONVERT 

TO 

INTEGER 

NAVDA  = 

AVDA  +  .49 



COMPUTE   NO.  OF  DAYS  BETWEEN 
IDATE   AND  MATURITY  DATE 
JDIF  =  J  DATE  +  NAVDA-  I  DATE 
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Figure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 

® 



COMPUTE   STANDARD   DEVIATION  TO 
USE  IN  NORMAL  DISTRIBUTION  CURVE 
SIGMA  =  (  AVDA  /50. )  *  5.47 


COMPUTE   PLUS  NORMALIZED  Z 
VARIABLE   FOR  UPPER  INTEGRATION 
LIMIT 

ZP  =  ABS((FL0AT(JDIF))  +  .5)/SIGMA)/V2" 


TEST  FOR 
INSIGNIFICANT  YIELD  CONTRlBUTIONT 
IZPI    >  2.7  ? 


YES 


NO 


COMPUTE  MINUS  NORMALIZED  Z  VARIABLE 
FOR  LOWER  INTEGRATION  LIMIT      ZM  = 

ABS((  FLOAT  ( J  DIF  )) -.5 )/ SIG M  A  )  /  V2~" 


SCALE  ZP  TO  OBTAIN  SUBSCRIPT  FOR 
LOCATING  ERROR  FUNCTION  VALUE  IN 
TABLE 




LINEARLY  INTERPOLATE  ERROR  FUNCTION 
VALUE 


SCALE   ZM  TO  OBTAIN  LOWER  INTEGRATION 
SUBSCRIPT  AND  INTERPOLATE  IN  ERROR 
FUNCTION  TABLE 


0 
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Figure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


COMPUTE  AREA  UNDER  NORMAL  PROBABILTY 
CURVE  WHICH  IS  THE  FRACTION ,  F,  OF  BOLLS 
SET  ON  JDATE   MATURING  ON  IDATE 


HAVE  INTEGRATION   FROM  0  TO  .5  DAYS 
SO    F  =  EZP 

N 

^  

ACCUMULATE  CONTRIBUTION  FROM  JDATE 
TO  THE  YIELD  ON  IDATE 

YIELD  =  YIELD+  BL  SET  ( J  DATE  ) F 
C  YIELD  =  C  YIELD +  CBL  SET  IJ  DATE  )*  F 


COMPUTE  YOUNGEST  AND  OLDEST  AGES  OF 
SUSCEPTIBILITY  OF  THE  GREEN  BOLLS  TO 
PINK  BOLLWORM 

lYNG  =  10.*  (  AVDA /45.  )*.49 
lOLD  =  28.*(AVDA/45. )  ■*-.49 


IGRN   =  lYNG 


COMPUTE  AGE  OF  BOLLS  ON  IDATE 
JDIF  =  J  DATE  +  IGRN  -  I  DATE 


COMPUTE  STANDARD  DEVIATION  TO  USE  IN 
NORMAL  DISTRIBUTION  CURVE 
SIGMA  =  (FLOAT  {IGRN  )/50.  5.47 
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Figure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


© 


COMPUTE  PLUS  NORMALIZED  Z  VARIABLE 
ZP  =  ABS  ((FLOAKJDIF  +  .  5)  /  SIGMA)/ V2~ 


TEST  FOR 
INSIGNIFICANT  YIELD  CONTRIBUTIOI 
|ZP|  >  2.7  ? 


YES 


COMPUTE  MINUS  NO  R  M  A  LI  Z  E  D  Z  VA  R  I  A  B  LE 

FOR  LOWER  INTEGRATION  LIMIT 

ZM=  ABS({FL0AT{JDIF))-.5  )/SIGMA/V2~ 


SCALE  ZP  TO  FIND  SUBSCRIPT  AND 
INTERPOLATE  IN  ERROR  FUNCTION  TABLE 


SCALE  ZM  TO  OBTAIN  LOWER  INTEGRATION 
SUBSCRIPT  AND  INTERPOLATE  IN  ERROR 
FUNCTION  TABLE 


 ±  

COMPUTE  ERROR  UNDER  NORMAL  PROBABILTY 
CURVE  WHICH  IS  FR  ACT  I  ON  ,  F,  OF  BOLLS 
STARTED  ON  JDATE  WHICH  ARE  IGRN 
PHYSIOLOGICAL  DAYS  OLD  ON  IDATE 


Figure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


NO 


F  =  EZ  P 


ACCUMULATE  CONTRIBUTION  FROM  J  DATE 
TO  THE  NUMBER  OF  I6RN  DAYS  OLD  BOLLS 
PRESENT  ON  I  DATE 

GBLS(IGRN)  =  GBLS(I6RN)+  BLSET(JDATE)*F 
C6BLS(IGRN)=CBLS(IGRN)  +  CBLSET{JDATE)*F 


I6RN  =  IGRN+  1 


NO 


■^10 


Yl 

V 

:s 

JDATE  =  JDATE  +  1 

JDATE  >  IDATE  ? 


NO 


ACCUMULATE  HARVEST  OF  MATURE  BOLLS 

HRVST  =  HRVST  +  YIELD 
CHRVST  =  C HRVST +  CYIELD 
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Figure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


GET  INFESTATION  OF  PBW  IN  TERMS  OF  RATIO 
OF  NO. OF  LARVAE  TO  NO.  OF  GREEN  BOLLS 
WHEN  NO  CHEMICAL  TREATMENT  IS  USED 
CALL   INFEST  (IDATE ,  FNFST ) 


NO 


YES 


GET  AMOUNT  OF  NEW  INFESTATION  IN  YOUNGEST 
BOLLS  AS  IF  THEY  HAD  AVERAGE 
PNEW=FNFST:^GBLS(IYNG) 


NO 


PNEW  =  0 


NO 


YES 


WRITE   ERROR  MESSAGE 


IPl  =  1YN6  +  1 
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ure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


 ±  

IGRN  =  IPl 

^  

 ^  

SUM  THE  AMOUNT  NEW  INFESTATION  REPRESENTED 
BY  THE  CHANGE  IN  AVERAGE  INFESTATION  IN  THE 
GREEN  BOLLS  THAT  HAVE  BEEN  SUSCEPTIBLE 
MORE  THAN  ONE  DAY. 

PNEW  =  PNEW  +  (FNFST-FYEST  )^GBLS  (IGRN) 


 ±  

IGRN  =  IGRN  +  1 


 T 

UPDATE   YESTERDAY  VALUES 
FYEST  =  FNFST 
lYYEST  =  lYNG 


TGBLS 

=  0 

TCGBLS 

=  0 

IGRN 

=  lYNG 

 ±  

SUM  TO  GET  THE  TOTAL  NO.  OF  GREEN  BOLLS 
W/0  AND  W/  CHEM  TMT. 

T6BLS=  TGBLS  +  GBLS(IGRN) 
TCGBLS =  TCGBLS 4  CGBLS(IGRN) 


 ±  

IGRN  =  IGRN  +  1 


NO 
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gure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield 
pink  bollworm  model. — Continued 


INITIALIZE  AMOUNT  OF  NEW  INFESTATION 
WITH  CHEMICAL  TREATMENT 

CPNEW  =  0 


 ^  

COMPUTE  AVERAGE  NUMBER  OF  NEW 
LARVAE  PER  GREEN  BOLL  WHEN  NO 
CHEMICAL  TREATMENT 

PORTN  =  PNEW/TGBLS 


IGRN  =  lYNG 

^  

 i  

COMPUTE  THE  NO.  OF  PBW'S  PER  BOLL  PER 
PARTICULAR  DAY'S  AGE  FOR  NO  CHEM  TMT 
ASSUMING  THAT  THESE  NEW  LARVAE  INFEST  ALL 
SUSCEPTIBLE  BOLLS  WITH  EQUAL  PREFERENCE 
PVBVD  =  PORTN*GBLS  (IGRN  ) 

 ±  

COMPUTE  THE  AMOUNTOFNEW  INFESTATION 
EXPECTED  WHEN  THERE  IS  A  CHEMICAL  TMT. 

CPNEW  =  CPNEW+PVBVD-3j(-CGBLS(IGRN)/GBLS 

(IGRN) 


 ±  

IGRN  =  IGRN  +  1 


GET  FRACTION  OF  LARVAE  GOING  INTO 
DIAPAUSE  FROM  TODAY'S  INFESTATION 
CALL  DIPAS  (IDATE,  FDIPAS  ) 
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Figure  12. — Flow  diagram  for  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


ACCUMULATE  NUMBER  OF  LARVAE  GOING  INTO 
DIAPAUSE  W/0  AND  W/  CHEM  TMT 

DPBW  =  DPBW  +  PNEW*FDIPAS 
CDPBW  =  CDPBW  +  CPNEW-*  FDIPAS 


OUTPUT 

DAY,  MONTH  ,  IDATE, MATURITY  DATE, 
FNFST,  FDIPAS,  BLSET  (IDATE),  CBLSET 
(IDATE  )  ,  T6BLS,TC6BLS,HRVST,CHRVST, 

DPBW , CDPBW 


IDATE  =  IDATE+  1 

COMPUTE  HARVEST  AND  PBW  %  REDUCTION 
DUE  TO   CHEM  TMT. 

HRDUCT  =  (HRVST-  CHRVST  )  *  100.  /  HRVST 
PRDUCTMDPBW-CDPBW-H^IOO./DPBW 


OUTPUT 

HRDUCT 
PR  DUCT 


s 

IRUN  =  IRUN  +  1 

NO 
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Figure  13a. — Computer  listing  of  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. 


CCri  rON   YIELD   -   PINK   BOLLWORM  MODEL 


PROGRAM  TO  COMPUTE  THE  COTTON  YIELD  AND  DIAPAUSINO 
PINK  BOLLWORM  POPULATIONS  FOR  FIELDS  THAT  HAVE  AND  HAVE  NOT 
BEFN  TREATED  WITH  CHEMICALS  TO    TERMINATE  LATE  SEASON  BOLL  SE 


PROGRAMMED  BY  B.     A     KIMBALL    (MAY,  1976) 
U.    S.    WATER  CONSERVATION  LABORATORY 
4331   E.  BROADWAY 
PHOENIX  AZS5040 


FOR  DATA  GENERAL  CORPORATION  NOVA  SCO  JUMBO 
MINICOMPUTER  WITH  32K  OF  MEMORY,    CRT   TERMINAL,    LINE  PRINTER, 
AND  DISK. 


D I  MENS  I  ON  BLSE  l'  ( "iiOO  )  ,  CBLSE  l'  (  500  )  ,  GBLS  (  1 00  )  ,  CGBLS  (  1 00  ) 
D I MENS  I UN   I T 1 TLE ( 40 ) ,  ND AT  E ( 3 ) , E ( 2S0 ) , KCSTOR ( 20 ) 
TYPE   "  COTTON  YIELD  -  PINK  BOLLWORM  MODEL" 

ACCEPT  "OUTPUT  TO  LPT  (12)    OR  C0TRSLT(2)  ",IOUT 
IF( lOUT.  NE.  12)    CALL  OPEN   ( lOUT, "COTRSLT" ) 
TYPE   "RUN   INFORMATION  V  " 
READ (11,2)    I T I TLE ( 1 ) 
FORMAT ( S72 ) 

ACCEPT  "NO.    OF   RUNS   ",  NRUNS 
DO  23   IRUN=1 ,  NRUNS 

ACCEPT""CHEMICAL  TREATMENi'   DATE  (D,  M)  ",ID,MO 

CALL  DIES(  ID,  MO,  75) 
KCSTOR ( IRUN)=ID 
CONTINUE 

FIRST   FLOWERING  DATE 
KFLWR=154 

FINAL  HARVEST  DATE 
KHRVST^334 


START  RUN  LOOP 

DO  9999  IRUN=1,NRUNS 

KCHEM=KCSTOR ( I RUN ) 


WRITE   RUN  INFORMATION 
WR I TE ( I  OUT ,15)    I  T  I T  LE ( 1 ) 
FORMAT ( IX, S72) 
WR I TE ( I OUT , 16)  KCHEM 

FORMAT  ("   DATE  OF  CHEMICAL  TREATMENT  =   ",I4./"  ") 
GFT  AND  WRITE  RUN  TIME 
CALL  FGTlMdHR,  IMIN,  I  SEC) 
CALL  DATE(NDATE,  lER) 

WR IT  E ( 1 OUT  ,  13)    NDATE ( 2 )  ,  ND ATE ( 1 )  ,  ND ATE ( 3 )  ,  I HR ,  IMIN,  I SEC 
FORMAT  (  "    RUN  ON    ",3I3,  "      AT    ",3I3./"  ") 
WRITE  HEADINGS 
WRIT  E( luUT  ,  14) 

FORMAT (zSX,  "BOLL  SET     GREEN  BOLLS         HARVEST     DIAPAUSE  PBW". 
1  2  X  ,  "  M  A  TR  "  ,  IX,  "I  F  S  D I  PS  "  ,  4  (  "        CK       T  MT    "  ) ./ 
"   DA  MO  DATE  DATE   ",2("   7.      "  )  ,  5X ,  " ./A/DA "  , 
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Figure  13b. — Computer  listing  of  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


3  4X,  8(  "        )  ,  "/.  01  A",  Si  "   -"  )  ) 

SET  UP  ERROR  FUNCTION  TABLF 
[10  30   1  =  1,  2S0 
2=FL0AT(I-1 )/100 
E( 1 )=ERFUNv  Z ) 

0  CONTINUE 

INITIALIZE  HARVEST  AND  NUMBER  OF  PINF  BOLLWORMS   IN  DIAPAUSE 

HRVST=0. 

CHRVST=0. 

DPBW=0. 

CDPt:W=0. 

INITIALIZE  YESTERDAY'S  YOUNGEST  SUSEPTIDLE  ACE  AND  YESTERDAY ■- S 

AVERAGE  FRACTIONAL  INFESTATION 

IYYEST=9 

FYES7-0 

INITIALIZE  DOLLSET  ARRAYS 
DO   1    1=1, KHRVST 
ELSE  I  ( I ) =0 
CBLSE  I  (  I  )  =0. 

1  CONTINUE 

COMPU  I  E  BOLL  SE  [   WITHOU  I"  AND  WITH  CHEMICAL  TREATMENT, 

DO  3   I DATE=KFLWR, KHRVST 

CALL   BOLSE  r ( I DATE ,  B ) 

BLSET( IDATE)=B 

JDIF=IDATE-KCHEM 

CALL  CHEM ( JD I F .  FR AC  1  ) 

CBLSE  r  (  I  DA  rE  )  =B-«-FRAC  T 

3  CONTINUE 

STARl    OUTER  LOOP 

DO   999    I DATE=KFLWR, KHRVST 

INITIALIZE  YIELDS  OBTAINED  ON   I DATE 

YIELD=0 

CYIELD=0. 

INITIALIZE   ARRAYS   OF  GREEN   BOLLS   PRESENT   ON    I  DATE 
DO  4   1=1, 100 
GBLS( I )=0 
CGBLS  (  I  )  =0. 

4  CONTINUE 

START    SECONDARY   LOOP   TO  CONTRIBUTION  FROM  ALL   PREVIOUS   DATES  YIELD 

AND  GREEN   BOLL  ARRAY 

DO  998  .JDATE=KFLWR,  I  DATE 

GET   NO     OF  DAYS  THE  BOLLS  SE V  ON  JDATt  WILL    1 AKE  TO  MATURE 
CALL   BOLMAK  JDATE,  AVDA  ) 
NAVDA=AVDA+.  49 
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Figure  13c. — Computer  listing  of  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 

GET  CONTRIBUTION   FROM   JDATE  TO   THE   YIELD   ON    I DATE 
COMPUTE   NO     OF   DAYS  FROM   MATURITY   DATE   10    I  DATE 
JD I F  =  JDAT  E  +  NAVDA- I  DATE 

COMPUTE  STANDARD  DEVIATION  TO  USE   IN  NORMAL  DISTRIBUTION 
CURVE   AND   THE   NORMALIZED  VARIABLE.  2. 
S  JGMA=<  AVDA/50.  ) -fi-S.  47 

ZP= (  ( ABS ( FLOAT ( JD I F )  ) +   5 ) /S I GMA ) / 1 .  41 42 1 
I F ( ABS ( ZP ) .  GT.  2  7 )    GOTO  20 

Z M= (  ( ABS ( FLOAT ( JD IF))-.  5 ) /S I GMA ) / 1 .  41421 

GFT   FRACTION  OF  THE   BOLLS  STARTED  ON   JDATE   MATURING   ON    I  DATE 

XP^-ZP*100.       +  1. 

iP=XF 

EZP=E<  IP)  +  (  XP-FLOAT(  I P  )  )    (  E  (  I P+ 1  ) -E  (  I P  )  ) 
IF(JDIF.  EO.  0)    GOrO  22 
XM=ZM*100.  +1 
IM=XM 

EZM=E( IM)+( XM-FLOAT( 1 M ) ) * ( E ( I M+ 1 ) -E ( I M ) ) 

F=(EZP-EZM)/2. 

IF(JDIF.  EQ.  0)  F=E2P 

Y I  ELD= Y I ELD+ELSET  (  JDATE  )  -s-F 

CY  I  ELD=CY  I ELD+CBLSE T  (  JDATE  )  -s-F 


GE'I'  CONTRIBUTION   FROM   JDATE  TO   THE   NUMBER  OF  GREEN  BOLLS  FOR 
EACH   DAY   OF   AGE   THAT   THE   BOLLS   ARE  SUSEP'TIELE   TO   PINK  BOLLWORM 
FIRST   NORMALIZE  FIRST    AND  OLDEST  SUSEPT I B I L I TY   AGES  FOR  THE 
CHANGE   IN  MATURATION  AGE  DURING  SEASON 
IYNG=10.  ■fi-(AVDA/45    ) +.  49 
I0LD=28.       AVDA/4^;l.  )    +  .  49 
no  997   IGRN=IYNG, I OLD 

GET  CONTRIBUTION   FROM   JDA'TE  TO  THE   NUMBER   OF    I  CRN   DAYS   OLD  BOLLS 

PRESENT   ON    I  DATE 

JD I F= JDA  rE+ I GRN- I  DATE 

COMPUTE  STANDARD   DEVIATION   TO   USE    IN  NORMAL   DISTRIBUTION  CURVE 

ANTi  NORMALIZED   VARIABLE.  Z 

S I  GMA=  (  FLOAT  (  I  GRN  )  /50.  )  47 

2P= (  ( ABS ( FLOAT ( JD I F ) ) +.  5 ) /S I GMA ) / 1 .  41 42 1 

IF(  ABS(  ZP).  GT.     2  7)    GOTO  997 

7M=  V  (  ABS  (  FLOAT  (  JD  I F  )  )  -.  5  )  /S I  GMA  ) ./ 1 .  41421 

GFl    FRACTION  OF   "THE   BOLLS   STARTED  ON   JDATE  WHICH     ARE    I GRN 
DAYS  OLD  ON  IDATE 
XP=ZP-ft-100.       +  1. 
IP=XP 

EZP=E ( i  P )  +  ( XP-FLOAT ( I P )  ) * ( E ( I P+ 1 ) -E ( I P ) ) 
IF(JDIF.  EQ.  0)    GOTO  21 
XM=ZM-&100.     +  1. 
IM=XM 

E  Z  M=E ( I M )  +  ( X  M-FLO AT (  I M )  ) * ( E ( I M+ 1 ) -E ( 1 M )  ) 
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Figure  13d. — Computer  listing  of  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. --Continued 


F-(I^ZP-EZM)/2. 

II-  (  JDIF  EQ.  0)  F=EZP 

GBLS  (  I GRN  )  =GBLS  (  1 GRN  )  +ELSET  (  JDA'I  E  )  *F 
CGBLS  (  1 HRN  )  =CGBLS  (  I GKN  )  +i::BLSET  (  JUATE  )  *F 

CONTINUE 
CONTINUE 

ACCUMULATE  HARVEST  OF   MATURE  BOLLS 
HRVST=HRVST+Y1ELD 
CHRVST=CHRVST+CY I ELD 

COMPUIE  THE  AVERAGE  RATIO  OF  LARVAE  PER   "10-2S"   DAY  OLD  BOLLS  UNDER 
NORMAL  CONDFflONS  WHEN  NO  CHEMICAL    rREATNENI    IS  USED 
CAI  .L    1NPEST(  I  DATE,  FNFST) 

GET  THE  AMOUT  OF  NEW   INFESTATION  FOR  NO  CHEN  TMT. 
I F ( 1 YNG.  EQ.  I YYEST )    PNEW=FNFST#GBLS ( I YNG ) 
IF(1YNG  EQ   IYYEST+1)  PNEW=0. 

IFdYNG  NE   lYYESr  AND   lYNG  NE   IYYEST+1)    WR I TE  v  I  OUT,  5 )  IDATE 

FORMAT ( "   GREEN  BOLLS   INITIAL  DAY  OF  SUSEPT I B I L I TY  ALIGNMENT  ERROR"/ 

"    IDATE  =  15) 

IP1=IYNG+1 

FDIF=FNFST-FYEST 

DO  A   IGRN=IP1, I OLD 

PNEW=PNFW+FD  1  F*-GBLS  (  I  GRN  ) 

CONT I NUE 

UPDATE  YESTERDAY  VALUES 

FYEST=FNFST 

IYYEST=IYNG 

GFl    THE    rOTAL  NUMBER  OF  GREEN  BOLLS 

TGBLS=0, 

TCGBLS=0. 

DO  7   1GRN=1YNG, I OLD 
TGELS=  i  GBLS+GBLS ( I GRN ) 
TCGBLS=TCGBLS+CGBLS( IGRN) 
Ci:  IN  T  I  NUE 

STAR!    LOOP  FOR  C0MPU1ING  THE   INFESTATION  WHEN  THERE   IS  A  CHEM  TMT" 
CPNEW=0. 

PORTN=PNEW/ IGBLS 
DO  S   1GRN=IYNG, I OLD 

COMPUTE  THE  NUMBER  OF  LARVAE  PER  BOLL  PER  PARTICULAR  DAY  S  AGE  FOR 
NO  CHEM  TMT',    ASSUMING  THAT   THESE  NEW  LARVAE   INFEST'  ALL  SUSEPT  I BLE 
BOLi S  WITH  EQUAL  PREFERENCE. 
PVBVD=PORTN-«-GBLS  (  IGRN ) 

NOW  KNOWING  THE  LARVAE/BOLL/DAY  AGE  FOR  THE  IGRN  1 H  DAY  AGE  FOR  NO 
CHFM  TMT,  COMPUTE  THE  AMOUNT  OF  NEW  INFESTAION  EXPECTED  WHEN  THERE 
IS  A  CHEM  TMT. 

CPNEW=CPNEW+PVBVD^<-CGBLS(  1  GRN). /GBLS  (  IGRN) 
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Figure  13e. — Computer  listing  of  the  main  program  of  the  cotton  yield- 
pink  bollworm  model. — Continued 


S  CONTINUE 

CUMPUTE  FRACTION  OF  LARVAE  GOING   INTO  DIAPAUSE  FROM  TODAY -"S 
INFE0I  A  riON. 

CALL   D I PAS ( I DATE. FD I PAS ) 

ACCUMULATE   NUMBER   OF   LARVAE   GOING    I N 1 0   DIAPAUSE  W/0   AND  W/  CHEM  TMT. 
DPBW=DPBW  +  PNEW-K-FDIPAS 
CDPBW=CDPBW+CPNEW*FD I PAS 

OUTPUT  RESULTS 
IO=NAVDA+IDATE 
I  l=FNFS  l'-ft-100.         +.  49 
I2=FDIPAS-frl00.         +  49 
I3=ELSE1( I  DATE)    +.  49 
I4=CBLSE  I'(  IDATE)    +  49 
15=1 GBLS/ 100     +  49 
I6=TCGBLS/100.     +.  49 
I7=HRVST/100.     +.  49 
IS=CHRVST/100.     +.  49 
I9=DPBW/100.     +.  49 
1 10=CDPEW/ 100.     +.  49 
IEi=IDATE 
M0=0 

CALL  DIES(  ID,  MO,  Tlf.) 

WR1TE<  lOUT,  11  )  ID,  MO,  IDATE,  10,  II,  12,  IS,  14,  15,  16,  17,  IS,  Iv,  HO 
1  1        FORMAT  (  IX,  2(  12,  ",  "  )  ,  2(  14,  ",  "  )  ,  2(  IS,  "  ,  "  )  ,  7  (15,  ",  "  )  ,  15) 

999  CONTINUE 

COMPUIE  HARVEST  AND  DIAPAUSED  PINK 
CHEMICAL  TREATMENT. 

HRDUC  r=--  (  HRVST-CHRVST  )  *  1 00.  ./'HRVST 
PRDUCT=  (  DPBW-CDPBW  )  *  1 00.  ./DPBW 
WR1TE(  lUUl",  12)    HRDUCT,  PRDUCT 
12        FORMAl  (  "      "./■"    HARVEST  REDUCTION 
1  "   PINK  BOLLWORM  REDUCTION 

DO  24   1  =  1,  6 

24  WRITE ( I OUT, 25) 

25  FORMAT (  "  "  ) 
9999  CONTINUE 

STOP 
END 


Figure  13f. — CUBIC  function. 
Fl.  INCT  I  ON   l:UB  I C  (  A ,  B ,  C ,  D ,  X  ) 

CUBIC   EVALUATES  THE  CUBIC  POLYNOMIAL 
CUBIC=   A   +   B-s-X    +   C-«-X-s--it-2  +  D-!i-X«--«-3 
GIVEN   A,  B,  C,  D,    AND  X. 

CUBIC=   A  +   X-ft-(B   +    X-fi-v'C  +  X-R-D)) 

RFIURN 

END 


BOLLWORM  REDUCTION  DUE  TO 


(X)  =  ",  F6  1./ 
iX)    =   " ,  F6.  1  ) 
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Figure  13g. — ERFU1^1  function. 


COMPILER  DOUBLfc  PRECISION 
FUNCTION  ERFUN(X) 

FUNCTION  ERFUNvX)    COMPUTES:  THE  ERROR  FUNCTION  FOR  A  OIVEN 
VALUE  OF    X   FOLLOWING  THE  SERIES  590.    P   136   IN  UWIGHT,    B  H  , 
TABLES  OF    INTEGRALS   AND  OTHER  MATHEMATICAI     DATA,    MACNlLLAN  CO.  , 
N  Y.  ,  196L 

TEST  FOR  LARGE  X 
FRFUN=-1 

IF(  X,  LF.  -4.  )  RETURN 
FRFUN=+ 1 

IF ( X   GE  +4    ) RET URN 

I N IT  I AL  I  ZE  SUM,  S  1 GN,  NUMFRAT  OR,  F  ACT  OR  1 AL 

SUM=  1 

.JSIGN=1 

FACT  =1 

XSQ=X-»X 

XNUN=1. 

START   ITERATION  LOOP 

N=0 

N=N+1 

F  ACT  =F  ACT  ^^FLOAT  (  N  ) 

DENOM=FACT*FLOAT  (  2-P-N+  1  ) 

XNUM=XNUM*XSQ 

JSIGN=.JSIGN-R-(-l  ) 

T  ERN-F  LOAT  ( JS 1 GN ) * XNUM/ DENOM 

SUM=SUM+TERM 

IF(ABS(TERM)    LE   1   E-S )    GOTO  2 
GCTTO  1 

CONVERGENCE  AT  T  A  1 NED 
FRFUN=X^<-SUM-«- 1    1 28379 
RET  URN 
END 


Figure  13h. — CHEM  subroutine  for  2,4-D  plus  chlorf lurenol  (fig.  5). 

Z,  4-D  +  CHLORFLURENOL 
SUEROUT  I NE   CHEM    ( JD 1 F ,  F  PACT  ) 

SR  CHEM  COMPUTES   THE   FRACTIONAL   REDUCTION    IN  BOLL  SET, 
PPACT,    AS  A  FUNCTION  OF  DAYS  FROM  TREATMENT,    JDIF,  FOR 
COTTON  TREATED  Wl TH  2, 4-D  +  CHLORFLURENOL. 

X^JDIF 

IF (JDIF  LF    -10)  FRACT=100. 
IF  (  JDIF   G  1  .  -10     AND.     JDIF   LF.  23) 
1  FRACT=CUBIC(.  44'D75bE2,  -  39977SE1,  .  181024,  -  40':.049E-2,  X) 

IF  (JDIF.  GT.  23)  FRACT=0. 
FRACT=FRACT./  100. 
RET  URN 
END 
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Figure  13i. — CHEM  subroutine  for  TD-1123  plus  chlorf lurenol  (fig.  6). 


1111123  +  CHLORF  LURENOL 
SUBROUT I NE   CHEM  <  JD 1 F .  F  R ACT ) 

SR  CHEM  COMPUIES   THE  FRACl'IONAL   REDUCTION    IN  BOLL  SET 
,    FRACT,    AS   A  FUNCTION   OF   DAYS   FROM  TREATMENT,    JLHF,    FOR  COITON 
TREATED  WITH   TD1123   +  CHLORFLURENOL. 

X=JDIF 

IF(JD1F.  LE.  -10)  FRACr=100^ 
IFCJDIF.  GT   -10   .  AND     JDIF.  LE.  -5) 

FRAC1=CUBIC(-.  V47711E7.  -.  5647b4E2,  -.  109v052E2,  -.  76663^,  X) 
IF  (JDIF.  GT.  -5   .  AND.     JDIF.  LF  -3) 

FRACT=CUBIC( -.  155346E1,  -.  544S24,  .  280S92,  -.  208937E-1,  X) 
IF  (JDIF.  GT.  -3)  FRACT=0 
FRACl  =FRACT./100 
RE"(  URN 
END 


Figure  13j.~B0LSET  subroutine  for  DPI  16  cotton  (fig.   1) . 

5973  DEL     16  COTTON     RARKER  AZ 
SUBROUT  I  NE  BOLSET  (  1  DAT  E ,  BOLI  .S  ) 


BOLLS,  THE  NUMBER  OF  BOLLS 
riON   OF    I DATE,    THE  GIVEN 


THIS  SR  BOLSET  COMPUTE 
SET  PER  ACRE  PER  DAY,    AS  A  FUNC 
JULIAN  DATE   THROUGH   THE   GROWING  SEASON  FOR   DPL    16  COTTON 
IN    1973   AT"   PARKER,  AZ. 


X=IDATE 

IF (I  DA  IE  LE   163)  BOLLS=0. 

IFdDATE.  GT.  163  .  AND.     I  DATE.  LE.  175) 

BOLLS=CUBIC(-.  219660ES,    392S13E6,  -.  234316E4,    46626SE1,  X) 
IF(IDATE.  or.  175   .  AND.     I  DATE   .  LE.  1S3) 

BOLLS=CUBIC(.  736304ES,  -.  124593E7,  .  702139E4,  -.  131746E2,  X) 
IFdDAT  E.  GT'.  1S3  .  AND.     I  DATE.  LF.  187) 

B0LLS=CUB1C(-.  112125E9,.  1 799 1 9E7, -.  96 1 S90E4,  .  171356E2,  X) 
IFdDATE  or.  187   .  AND.     I  DATE.  LF.  195) 

BOLLS=CUBICv-.  515414E7,  .  S30775E5,  -.  441827E3,  .  777208,  X) 
IFdDATE  GT.  195  .  AND      I  DATE.  LE  288) 

B0LLS=CUB1C(.  121532E7,  -.  149142E5,  .  606948E2,  -.  S18026E-1,  X) 
I F  (1  DATE  o  r.  288  )  BOLLS=0. 
RET  URN 
END 
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Figure  13k. — EOLSET  subroutine  for  AZ  64  cotton  (fig.  2). 


197":;  AZ64  LOTTuN  CRC  PHOENIX 

SUtiRUU'f  INt  BOLSt  f  (  IDA1E,.  BOLLS) 

THIS  SR  BOLStT   COMRUfLS  BOLLS,    THE  NUMBER  OF   BOLLS  SEV  PER 
ACRE   PER   DAY,    AS   A  FUNCTION   OP    I  DATE,    THE   OIVEN   JULIAN  DATE 
IHROUGH    THE   GROWING   SEASON   FOR  AZ64   COTTON    IN    1^'75  AT 
PHOENIX,  AZ. 


X^iDATE 

IF(lDArF.  LE    165)    BOLLS  =0. 

TF(1DA1E  GT.  165     AND.     IDATE  LF  190) 

B0LLS=CUB1C( .  670213E5-  -  10691 1E4 
IF<1DA'IE  GT.  190  .  AND      IDATE  LE  205) 

F0LLS=CUB1C(  27S513E6, -  440S53E4 
iP(  IDATE  GT.  205  .  AND.     IDAIE.  LF  215) 

B0LLS=CUB1C(-.  50506SE7,  .  73579SE5 
IFdDAlE  Gl.  215     AND      IDATE  LE.  230) 

BOLLS=CUBIC( .  S77021E7,  -  119270E6 
IF (IDATE  GT  230  .  AND      IDATE.  LF.  240) 

BOLLS=CUBIC( -.  769913E7,  955476E5 
IF  (IDATE.  GT.  240     AND.     IDATE  LF  260) 

BOLLS=CUBIC(-  226S2SE7,  276620E5 
IF (IDATE  GT  260  .  AND      IDATE  LF  300) 

B0LLS=CUB1C(    12526SE7, -.  129645E5 
IF  (  IDATE.  GT.  300)    BOLLS  =0. 
RP I  URN 
END 


551010E1, -   904104E-2, X) 

230S60E2, -.  39S760E-1,  X) 
-,  357  345E3,  .  57S7  1  1 ,  X  ) 

539630E3,  -.  81  194S,  X  ) 
-.  39435SE3,  .  541659,  X  ) 
-.  1  1  1502F3,  .  143302,  X  ) 

447539E2, -.  515275E-1, X ) 


Figure  131. — INFEST  subroutine  for  no  insecticide  (fig.  7). 

C:  NO   INSECTICIDE  TREATMENT 

SUBROU I  I NE   I NFEST ( I DAT  E ,  FNFS T ) 

C 

C  THIS  SR   INF ESI'  COMPUTES  THE  FRACTION  OF  GREEN  BOLLS  INFESTED 

C  WIT  H  PINK  BOLLWORN  LARVAE,    FNFS  I',    AS  A  FUNCTION  OF   IDATE,  THE 

C  GIVEN  JULIAN  DATE  THROUGH    IHE  GROWING  SEASON 

I": 

X^ I DATE 

IF ( IDATE  LE   196)  FNFST=0 

I F( IDATE  GT.  196  .  AND.     IDATE  LF  260) 

1  FNFST=CUBIC(-  797514E3,  .  139S97E2,  -.  S12771E-1,  .  1564S1E-3,  X) 
IF ( IDATE  GT  260     AND      IDATE.  LE  270) 

2  FNFST=CUB1C(    S29357E5, -  952163E3,  .  363469E1,  -   460753E-2,  X) 
IF  (IDATE  or.  270  .AND.     IDATE  LE.  290) 

3  FNFST=CUBIC(-  271272E5,  .  27075SE3,  -.  394644,  .  9S4191E-3,  X) 
IFdDAlE.  GT.  290)  FNFST=157 

FNFST=FNFST/100 

RETURN 

END 
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Figure  13m. — INFEST  subroutine  for  insecticide  (fig.  8). 


I NSECT I C  i  DE  TREAT  El  l 

SUBRuUf I NE    1 NFES1  ( 1  DATE,  F  NF  S7 ) 

"IHIS  SR  INFEST  COMPUTES  fHE  FRACTION  OF  GREEN  BOLLS 
INFESTEM  WITH  PINK  BOLLWORM  LARVAE.  FNFST;  AS  A  FUNCTION  OF 
11  (ATE,  THE  GIVEN  JULIAN  DATE  THROUGH  THF  GROWING  SEASON  FOR 
A  FIELD   THAT   HAS  EEFN  TREATED  WITH  INSECTICIDES. 

X=IDATE 

IF (I  DATE  LE.  210)  FNFST=0 

IF  (1  DATE.  GT.  210  .  AND.     iDATE  LE.  260) 

FNFST=CUBIC(-.  864694E3,  .  142253E-?',  -.  766563E-1,  .  135V04E-3,  X) 
IF  (  IDATE  GT.  260  .AND.     IDATE  LE.  270) 

FNFST=CUBIC( .  301939Eb,  -.  344143E3,  .  13016SE1,  -   163120E-2,  X) 
IF(  IDATE.  GT.  270) 

FNFST=CUBIC(-.  36S443E4,  .  322S2SE2,  -  924S6SE-1,  .  8999S6E-4,  X) 
FNFST=FNFST./100. 
RFTURN 
END 


Figure  13n.— BOLMAT  subroutine  (fig.  3). 
SUBROUT  INE   BOLMAT  v  IDATE,  AVDA ) 

SR   BOLMAT    COMPUTES   THE  NUMBER  OF   DAYS  FROM  FLOWERING 
TO   WHEN  50'/.   OF   THF   BOLLS   NATURF,    AVDA,    AS   A   FUNCTION   OF  THE 
JULIAN  DATE  OF   FLOWERING,  IDATE. 

X=IDATF 

AVDA=CUBIC(.  207076E2,  .  334656,  -.  271254E-2,  .  733372E-5,  X) 

RE  r URN 

END 


Figure  13o. — DIPAS  subroutine  (fig  9). 
SUBROUT  I  NE   D 1  FAS  (  1  DAT  E ,  FD  i  P At; ) 

SR   DIPAS   COMPUTES   THE   FRACTION   OF    PINK   BOLLWORM  LARVAE 
ENTERING   DIAPAUSE,    FDIPAS,    AS   A  FUNCTION   OF    JULIAN   DATE  THROUGH 
THE  GROWING  SEASON,  IDATE. 

X=IDATE 

IF  <  IDAT  E  LE.  235)  FDIPAS=0. 

IFv  IDATE.  or.  235  .  AND.     IDATE  LF  240) 

FDIPAS=CUBIC<-.  109131E6,  .  142503E4,  -.  620037E1,  .  39S959E-2,  X) 
IF  <  IDAT  E.  GT.  240     AND.     IDATE  LE.  260  ) 

FDIPAS=CUBIC(    1452S5E6,  -.  175517E4,  .  705045E1,  -  941433E-2,  X) 
IPC  IDATE   GT.  260      AND.     IDATE   LE   265 ) 

FDIPAS=CUBIC(-.  357154E6,    404221E4,  -.  152472F2,  .  191724E-1,  X) 
IF  (IDATE.  GT.  265  ,  AND.     IDATE  LE  300) 

FDIPAS=CUBIC<-.  430473E3,  .  3S2375E1,  -.  797S12E-2,  .  356943E-5,  X) 
I F  (  I  DAT  E.  GT.  300  )    F  D I P AS=95. 
FD  I  PAS=FD  I  PAS./"  1 00. 
RFTURN 
END 
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Figure  13p. — DIES  subroutine  for  calendar  conversion. 


SUBROU 1 1 NE  D 1 ES (ID, MO , I YR ) 

C 

C  CALENDAR  CONVERSION 

C  IE  MO  EQ  O,    THIS  SR  DIES  RETURNS  DATE  FROM  DAY  OF    YEAR  AND  YEAR 

C  IF  MO  NE  0,    THIS  SR  DIES  RETURNS  DAY  OF  YEAR  FROM  DAY,  MONTH, 

C  AND  YEAR,    LEAP  YEARS  CORRFCTED. 
C: 

C  FROM  G.    F     BYRNE,    AG     MET"  15(19/5)419. 


DIMENSION  MON( 12) 

MON( 1 )=31 

MON ( 2 ) =23 

MONO)  =31 

MON ( 4 ) =30 

M0N(5)=31 

MON ( 6 ) =30 

MON (7 )=31 

MON ( 3 ) =3 1 

MON ( 9 ) =30 

MON( 10) =31 

MON( 1 1 )=30 

MON( 12)=31 

IL=1 YR/4 

YL=FLOAT( lYR)/4 

I F  (  (  YL-FLOAT  (ID).  LT  .  .01)  MON  (  2  )  =MON  (  2  )  +  1 
I  =0 

IF  (  MO.  or.  O  )    GO  TO  4 

2  1=1+1 

I D= 1 D-MON ( I ) 
JF(ID-O)  3,3,2 

3  1D=1D+M0N(I) 
M0=  I 
RETURN 

4  1=1+1 
ID=ID+MON( 1 ) 

IF  (MO- I)  4,5,4 

5  :fD=ID-MON(l) 
M0=0 

RFT  URN 
END 
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